Isolation and characterization of a pangolin-borne HKU4-related coronavirus that potentially infects human-DPP4-transgenic mice.

We recently detected a HKU4-related coronavirus in subgenus Merbecovirus (named pangolin-CoV-HKU4-P251T) from a Malayan pangolin1. Here we report isolation and characterization of pangolin-CoV-HKU4-P251T, the genome sequence of which is closest to that of a coronavirus from the greater bamboo bat (Tylonycteris robustula) in Yunnan Province, China, with a 94.3% nucleotide identity. Pangolin-CoV-HKU4-P251T is able to infect human cell lines, and replicates more efficiently in cells that express human-dipeptidyl-peptidase-4 (hDPP4)-expressing and pangolin-DPP4-expressing cells than in bat-DPP4-expressing cells. After intranasal inoculation with pangolin-CoV-HKU4-P251, hDPP4-transgenic female mice are likely infected, showing persistent viral RNA copy numbers in the lungs. Progressive interstitial pneumonia developed in the infected mice, characterized by the accumulation of macrophages, and increase of antiviral cytokines, proinflammatory cytokines, and chemokines in lung tissues. These findings suggest that the pangolin-borne HKU4-related coronavirus has a potential for emerging as a human pathogen by using hDPP4.

Virus diversity, wildlife-domestic animal circulation and potential zoonotic viruses of small mammals, pangolins and zoo animals.

Wildlife is reservoir of emerging viruses. Here we identified 27 families of mammalian viruses from 1981 wild animals and 194 zoo animals collected from south China between 2015 and 2022, isolated and characterized the pathogenicity of eight viruses. Bats harbor high diversity of coronaviruses, picornaviruses and astroviruses, and a potentially novel genus of Bornaviridae. In addition to the reported SARSr-CoV-2 and HKU4-CoV-like viruses, picornavirus and respiroviruses also likely circulate between bats and pangolins. Pikas harbor a new clade of Embecovirus and a new genus of arenaviruses. Further, the potential cross-species transmission of RNA viruses (paramyxovirus and astrovirus) and DNA viruses (pseudorabies virus, porcine circovirus 2, porcine circovirus 3 and parvovirus) between wildlife and domestic animals was identified, complicating wildlife protection and the prevention and control of these diseases in domestic animals. This study provides a nuanced view of the frequency of host-jumping events, as well as assessments of zoonotic risk.

A SARS-CoV-2-Related Virus from Malayan Pangolin Causes Lung Infection without Severe Disease in Human ACE2-Transgenic Mice.

Coronavirus disease 2019 (COVID-19), which is caused by the novel coronavirus severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is the most severe emerging infectious disease in the current century. The discovery of SARS-CoV-2-related coronaviruses (SARSr-CoV-2) in bats and pangolins in South Asian countries indicates that SARS-CoV-2 likely originated from wildlife. To date, two SARSr-CoV-2 strains have been isolated from pangolins seized in Guangxi and Guangdong by the customs agency of China, respectively. However, it remains unclear whether these viruses cause disease in animal models and whether they pose a transmission risk to humans. In this study, we investigated the biological features of a SARSr-CoV-2 strain isolated from a smuggled Malayan pangolin (Manis javanica) captured by the Guangxi customs agency, termed MpCoV-GX, in terms of receptor usage, cell tropism, and pathogenicity in wild-type BALB/c mice, human angiotensin-converting enzyme 2 (ACE2)-transgenic mice, and human ACE2 knock-in mice. We found that MpCoV-GX can utilize ACE2 from humans, pangolins, civets, bats, pigs, and mice for cell entry and infect cell lines derived from humans, monkeys, bats, minks, and pigs. The virus could infect three mouse models but showed limited pathogenicity, with mild peribronchial and perivascular inflammatory cell infiltration observed in lungs. Our results suggest that this SARSr-CoV-2 virus from pangolins has the potential for interspecies infection, but its pathogenicity is mild in mice. Future surveillance among these wildlife hosts of SARSr-CoV-2 is needed to monitor variants that may have higher pathogenicity and higher spillover risk. IMPORTANCE SARS-CoV-2, which likely spilled over from wildlife, is the third highly pathogenic human coronavirus. Being highly transmissible, it is perpetuating a pandemic and continuously posing a severe threat to global public health. Several SARS-CoV-2-related coronaviruses (SARSr-CoV-2) in bats and pangolins have been identified since the SARS-CoV-2 outbreak. It is therefore important to assess their potential of crossing species barriers for better understanding of their risk of future emergence. In this work, we investigated the biological features and pathogenicity of a SARSr-CoV-2 strain isolated from a smuggled Malayan pangolin, named MpCoV-GX. We found that MpCoV-GX can utilize ACE2 from 7 species for cell entry and infect cell lines derived from a variety of mammalian species. MpCoV-GX can infect mice expressing human ACE2 without causing severe disease. These findings suggest the potential of cross-species transmission of MpCoV-GX, and highlight the need of further surveillance of SARSr-CoV-2 in pangolins and other potential animal hosts.

Induction of significant neutralizing antibodies against SARS-CoV-2 by a highly attenuated pangolin coronavirus variant with a 104nt deletion at the 3'-UTR.

SARS-CoV-2 related coronaviruses (SARS-CoV-2r) from Guangdong and Guangxi pangolins have been implicated in the emergence of SARS-CoV-2 and future pandemics. We previously reported the culture of a SARS-CoV-2r GX_P2V from Guangxi pangolins. Here we report the GX_P2V isolate rapidly adapted to Vero cells by acquiring two genomic mutations: an alanine to valine substitution in the nucleoprotein and a 104-nucleotide deletion in the hypervariable region (HVR) of the 3'-terminus untranslated region (3'-UTR). We further report the characterization of the GX_P2V variant (renamed GX_P2V(short_3UTR)) in in vitro and in vivo infection models. In cultured Vero, BGM and Calu-3 cells, GX_P2V(short_3UTR) had similar robust replication kinetics, and consistently produced minimum cell damage. GX_P2V(short_3UTR) infected golden hamsters and BALB/c mice but was highly attenuated. Golden hamsters infected intranasally had a short duration of productive infection in pulmonary, not extrapulmonary, tissues. These productive infections induced neutralizing antibodies against pseudoviruses of GX_P2V and SARS-CoV-2. Collectively, our data show that the GX_P2V(short_3UTR) is highly attenuated in in vitro and in vivo infection models. Attenuation of the variant is likely partially due to the 104-nt deletion in the HVR in the 3'-UTR. This study furthers our understanding of pangolin coronaviruses pathogenesis and provides novel insights for the design of live attenuated vaccines against SARS-CoV-2.

Exploring the Natural Origins of SARS-CoV-2 in the Light of Recombination.

The lack of an identifiable intermediate host species for the proximal animal ancestor of SARS-CoV-2, and the large geographical distance between Wuhan and where the closest evolutionary related coronaviruses circulating in horseshoe bats (members of the Sarbecovirus subgenus) have been identified, is fueling speculation on the natural origins of SARS-CoV-2. We performed a comprehensive phylogenetic study on SARS-CoV-2 and all the related bat and pangolin sarbecoviruses sampled so far. Determining the likely recombination events reveals a highly reticulate evolutionary history within this group of coronaviruses. Distribution of the inferred recombination events is nonrandom with evidence that Spike, the main target for humoral immunity, is beside a recombination hotspot likely driving antigenic shift events in the ancestry of bat sarbecoviruses. Coupled with the geographic ranges of their hosts and the sampling locations, across southern China, and into Southeast Asia, we confirm that horseshoe bats, Rhinolophus, are the likely reservoir species for the SARS-CoV-2 progenitor. By tracing the recombinant sequence patterns, we conclude that there has been relatively recent geographic movement and cocirculation of these viruses' ancestors, extending across their bat host ranges in China and Southeast Asia over the last 100 years. We confirm that a direct proximal ancestor to SARS-CoV-2 has not yet been sampled, since the closest known relatives collected in Yunnan shared a common ancestor with SARS-CoV-2 approximately 40 years ago. Our analysis highlights the need for dramatically more wildlife sampling to: 1) pinpoint the exact origins of SARS-CoV-2's animal progenitor, 2) the intermediate species that facilitated transmission from bats to humans (if there is one), and 3) survey the extent of the diversity in the related sarbecoviruses' phylogeny that present high risk for future spillovers.

The high diversity of SARS-CoV-2-related coronaviruses in pangolins alerts potential ecological risks.

Understanding the zoonotic origin and evolution history of SARS-CoV-2 will provide critical insights for alerting and preventing future outbreaks. A significant gap remains for the possible role of pangolins as a reservoir of SARS-CoV-2 related coronaviruses (SC2r-CoVs). Here, we screened SC2r-CoVs in 172 samples from 163 pangolin individuals of four species, and detected positive signals in muscles of four Manis javanica and, for the first time, one M. pentadactyla. Phylogeographic analysis of pangolin mitochondrial DNA traced their origins from Southeast Asia. Using in-solution hybridization capture sequencing, we assembled a partial pangolin SC2r-CoV (pangolin-CoV) genome sequence of 22 895 bp (MP20) from the M. pentadactyla sample. Phylogenetic analyses revealed MP20 was very closely related to pangolin-CoVs that were identified in M. javanica seized by Guangxi Customs. A genetic contribution of bat coronavirus to pangolin-CoVs via recombination was indicated. Our analysis revealed that the genetic diversity of pangolin-CoVs is substantially higher than previously anticipated. Given the potential infectivity of pangolin-CoVs, the high genetic diversity of pangolin-CoVs alerts the ecological risk of zoonotic evolution and transmission of pathogenic SC2r-CoVs.

Fundamental evolution of all Orthocoronavirinae including three deadly lineages descendent from Chiroptera-hosted coronaviruses: SARS-CoV, MERS-CoV and SARS-CoV-2.

The severe acute respiratory syndrome coronavirus (SARS-CoV) emerged in humans in 2002. Despite reports showing Chiroptera as the original animal reservoir of SARS-CoV, many argue that Carnivora-hosted viruses are the most likely origin. The emergence of the Middle East respiratory syndrome coronavirus (MERS-CoV) in 2012 also involves Chiroptera-hosted lineages. However, factors such as the lack of comprehensive phylogenies hamper our understanding of host shifts once MERS-CoV emerged in humans and Artiodactyla. Since 2019, the origin of SARS-CoV-2, causative agent of coronavirus disease 2019 (COVID-19), added to this episodic history of zoonotic transmission events. Here we introduce a phylogenetic analysis of 2006 unique and complete genomes of different lineages of Orthocoronavirinae. We used gene annotations to align orthologous sequences for total evidence analysis under the parsimony optimality criterion. Deltacoronavirus and Gammacoronavirus were set as outgroups to understand spillovers of Alphacoronavirus and Betacoronavirus among ten orders of animals. We corroborated that Chiroptera-hosted viruses are the sister group of SARS-CoV, SARS-CoV-2 and MERS-related viruses. Other zoonotic events were qualified and quantified to provide a comprehensive picture of the risk of coronavirus emergence among humans. Finally, we used a 250 SARS-CoV-2 genomes dataset to elucidate the phylogenetic relationship between SARS-CoV-2 and Chiroptera-hosted coronaviruses.

The hunt for SARS-CoV-2's ancestors heats up.

Cousins of the pandemic virus in animals help chart the path to humans.

SARS-CoV-2's claimed natural origin is undermined by issues with genome sequences of its relative strains: Coronavirus sequences RaTG13, MP789 and RmYN02 raise multiple questions to be critically addressed by the scientific community.

RaTG13, MP789, and RmYN02 are the strains closest to SARS-CoV-2, and their existence came to light only after the start of the pandemic. Their genomes have been used to support a natural origin of SARS-CoV-2 but after a close examination all of them exhibit several issues. We specifically address the presence in RmYN02 and closely related RacCSxxx strains of a claimed natural PAA/PVA amino acid insertion at the S1/S2 junction of their spike protein at the same position where the PRRA insertion in SARS-CoV-2 has created a polybasic furin cleavage site. We show that RmYN02/RacCSxxx instead of the claimed insertion carry a 6-nucleotide deletion in the region and that the 12-nucleotide insertion in SARS-CoV-2 remains unique among Sarbecoviruses. Also, our analysis of RaTG13 and RmYN02's metagenomic datasets found unexpected reads which could indicate possible contamination. Because of their importance to inferring SARS-CoV-2's origin, we call for a careful reevaluation of RaTG13, MP789 and RmYN02 sequencing records and assembly methods.

Molecular basis of cross-species ACE2 interactions with SARS-CoV-2-like viruses of pangolin origin.

Pangolins have been suggested as potential reservoir of zoonotic viruses, including SARS-CoV-2 causing the global COVID-19 outbreak. Here, we study the binding of two SARS-CoV-2-like viruses isolated from pangolins, GX/P2V/2017 and GD/1/2019, to human angiotensin-converting enzyme 2 (hACE2), the receptor of SARS-CoV-2. We find that the spike protein receptor-binding domain (RBD) of pangolin CoVs binds to hACE2 as efficiently as the SARS-CoV-2 RBD in vitro. Furthermore, incorporation of pangolin CoV RBDs allows entry of pseudotyped VSV particles into hACE2-expressing cells. A screen for binding of pangolin CoV RBDs to ACE2 orthologs from various species suggests a broader host range than that of SARS-CoV-2. Additionally, cryo-EM structures of GX/P2V/2017 and GD/1/2019 RBDs in complex with hACE2 show their molecular binding in modes similar to SARS-CoV-2 RBD. Introducing the Q498H substitution found in pangolin CoVs into the SARS-CoV-2 RBD expands its binding capacity to ACE2 homologs of mouse, rat, and European hedgehog. These findings suggest that these two pangolin CoVs may infect humans, highlighting the necessity of further surveillance of pangolin CoVs.

Exploration of hosts and transmission traits for SARS-CoV-2 based on the k-mer natural vector.

A severe respiratory pneumonia COVID-19 has raged all over the world, and a coronavirus named SARS-CoV-2 is blamed for this global pandemic. Despite intensive research into the origins of the COVID-19 pandemic, the evolutionary history of its agent SARS-CoV-2 remains unclear, which is vital to control the pandemic and prevent another round of outbreak. Coronaviruses are highly recombinogenic, which are not well handled with alignment-based method. In addition, deletions have been found in the genomes of several SARS-CoV-2, which cannot be resolved with current phylogenetic methods. Therefore, the k-mer natural vector is proposed to explore hosts and transmission traits for SARS-CoV-2 using strict phylogenetic reconstruction. SARS-CoV-2 clustering with bat-origin coronaviruses strongly suggests bats to be the natural reservoir of SARS-CoV-2. By building bat-to-human transmission route, pangolin is identified as an intermediate host, and civet is predicted as a possible candidate. We speculate that SARS-CoV-2 undergoes cross-species recombination between bat and pangolin coronaviruses. This study also demonstrates transmission mode and features of SARS-CoV-2 in the COVID-19 pandemic when it broke out early around the world.

Impacts of reference selection on the assembly of suspicious coronavirus genome.

The pandemic caused by SARS-CoV-2 has had a significant impact on the whole world. In a theory of the origin of SARS-CoV-2, pangolins are considered as a potential intermediate host. To assemble the genome of suspicious coronavirus (CoV) found in pangolins, SARS-CoV-2 was used as a reference in most of the previous studies, implicitly assuming the pangolin CoV and SARS-CoV-2 are the closest neighbors in evolution. However, this assumption may not be true. We investigated how the choice of reference genome affected the resulting CoV genome assembly. We explored various representative CoVs as the reference genome, and found significant differences in the resulting assemblies. The assembly obtained using RaTG13 as a reference showed better statistics in total length, N50, and pairwise distance reconstruction (PDR) scores than the assembly guided by SARS-CoV-2, indicating that RaTG13 may be a better reference. Therefore, RaTG13 should also be considered as a reference for assembling suspicious CoV found in pangolins and other potential intermediate hosts.

Urgent Need for Field Surveys of Coronaviruses in Southeast Asia to Understand the SARS-CoV-2 Phylogeny and Risk Assessment for Future Outbreaks.

Phylogenetic analysis of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is focused on a single isolate of bat coronaviruses (bat CoVs) which does not adequately represent genetically related coronaviruses (CoVs) [...].

Pathogenicity and transmissibility of a novel respirovirus isolated from a Malayan pangolin.

The identification of SARS-CoV-2-like viruses in Malayan pangolins (Manis javanica) has focused attention on these endangered animals and the viruses they carry. We successfully isolated a novel respirovirus from the lungs of a dead Malayan pangolin. Similar to murine respirovirus, the full-length genome of this novel virus was 15 384 nucleotides comprising six genes in the order 3'-(leader)-NP-P-M-F-HN-l-(trailer)-5'. Phylogenetic analysis revealed that this virus belongs to the genus Respirovirus and is most closely related to murine respirovirus. Notably, animal infection experiments indicated that the pangolin virus is highly pathogenic and transmissible in mice, with inoculated mice having variable clinical symptoms and a fatality rate of 70.37 %. The virus was found to replicate in most tissues with the exception of muscle and heart. Contact transmission of the virus was 100 % efficient, although the mice in the contact group displayed milder symptoms, with the virus mainly being detected in the trachea and lungs. The isolation of a novel respirovirus from the Malayan pangolin provides new insight into the evolution and distribution of this important group of viruses and again demonstrates the potential infectious disease threats faced by endangered pangolins.

Bat and pangolin coronavirus spike glycoprotein structures provide insights into SARS-CoV-2 evolution.

In recognizing the host cellular receptor and mediating fusion of virus and cell membranes, the spike (S) glycoprotein of coronaviruses is the most critical viral protein for cross-species transmission and infection. Here we determined the cryo-EM structures of the spikes from bat (RaTG13) and pangolin (PCoV_GX) coronaviruses, which are closely related to SARS-CoV-2. All three receptor-binding domains (RBDs) of these two spike trimers are in the "down" conformation, indicating they are more prone to adopt the receptor-binding inactive state. However, we found that the PCoV_GX, but not the RaTG13, spike is comparable to the SARS-CoV-2 spike in binding the human ACE2 receptor and supporting pseudovirus cell entry. We further identified critical residues in the RBD underlying different activities of the RaTG13 and PCoV_GX/SARS-CoV-2 spikes. These results collectively indicate that tight RBD-ACE2 binding and efficient RBD conformational sampling are required for the evolution of SARS-CoV-2 to gain highly efficient infection.

SARS-CoV-2 variants lacking ORF8 occurred in farmed mink and pangolin.

The SARS-CoV-2 Variant of Concern 202012/01 (VOC-202012/01) is rapidly spreading worldwide owing to its substantial transmission advantage. The variant has changes in critical sites of the spike protein with potential biological significance. Moreover, VOC-202012/01 has a mutation that inactivates the ORF8 protein, whose absence can change the clinical features of the infection. Why VOC-202012/01 is more transmissible remains unclear, but spike mutations and ORF8 inactivation stand out by their known phenotypic effects. Here I show that variants combining relevant spike mutations and the absence of ORF8 occurred in SARS-CoV-2 and related viruses circulating in other host species. A truncated ORF8 (Q23stop) occurred in a SARS-CoV-2-related virus from a pangolin seized in China in 2017, also with several mutations in critical spike sites. Strikingly, I found that variants without ORF8 (E19stop) and with the N501T spike mutation circulated in farmed mink and humans from Denmark. Although with differences to VOC-202012/01, the identification of these variants highlights the danger of having reservoirs of SARS-CoV-2 and related viruses where more transmissible variants may occur and spill over to humans.

Non-standard bioinformatics characterization of SARS-CoV-2.

A non-standard bioinformatics method, 4D-Dynamic Representation of DNA/RNA Sequences, aiming at an analysis of the information available in nucleotide databases, has been formulated. The sequences are represented by sets of "material points" in a 4D space - 4D-dynamic graphs. The graphs representing the sequences are treated as "rigid bodies" and characterized by values analogous to the ones used in the classical dynamics. As the graphical representations of the sequences, the projections of the graphs into 2D and 3D spaces are used. The method has been applied to an analysis of the complete genome sequences of the 2019 novel coronavirus. As a result, 2D and 3D classification maps are obtained. The coordinate axes in the maps correspond to the values derived from the exact formulas characterizing the graphs: the coordinates of the centers of mass and the 4D moments of inertia. The points in the maps represent sequences and their coordinates are used as the classifiers. The main result of this work has been derived from the 3D classification maps. The distribution of clusters of points which emerged in these maps, supports the hypothesis that SARS-CoV-2 may have originated in bat and in pangolin. Pilot calculations for Zika virus sequence data prove that the proposed approach is also applicable to a description of time evolution of genome sequences of viruses.

Evidence for SARS-CoV-2 related coronaviruses circulating in bats and pangolins in Southeast Asia.

Among the many questions unanswered for the COVID-19 pandemic are the origin of SARS-CoV-2 and the potential role of intermediate animal host(s) in the early animal-to-human transmission. The discovery of RaTG13 bat coronavirus in China suggested a high probability of a bat origin. Here we report molecular and serological evidence of SARS-CoV-2 related coronaviruses (SC2r-CoVs) actively circulating in bats in Southeast Asia. Whole genome sequences were obtained from five independent bats (Rhinolophus acuminatus) in a Thai cave yielding a single isolate (named RacCS203) which is most related to the RmYN02 isolate found in Rhinolophus malayanus in Yunnan, China. SARS-CoV-2 neutralizing antibodies were also detected in bats of the same colony and in a pangolin at a wildlife checkpoint in Southern Thailand. Antisera raised against the receptor binding domain (RBD) of RmYN02 was able to cross-neutralize SARS-CoV-2 despite the fact that the RBD of RacCS203 or RmYN02 failed to bind ACE2. Although the origin of the virus remains unresolved, our study extended the geographic distribution of genetically diverse SC2r-CoVs from Japan and China to Thailand over a 4800-km range. Cross-border surveillance is urgently needed to find the immediate progenitor virus of SARS-CoV-2.

Structure and binding properties of Pangolin-CoV spike glycoprotein inform the evolution of SARS-CoV-2.

Coronaviruses of bats and pangolins have been implicated in the origin and evolution of the pandemic SARS-CoV-2. We show that spikes from Guangdong Pangolin-CoVs, closely related to SARS-CoV-2, bind strongly to human and pangolin ACE2 receptors. We also report the cryo-EM structure of a Pangolin-CoV spike protein and show it adopts a fully-closed conformation and that, aside from the Receptor-Binding Domain, it resembles the spike of a bat coronavirus RaTG13 more than that of SARS-CoV-2.